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ABSTRACT 

Quasi-periodic, propagating fast mode magnetosonic waves in the corona were difficult to observe 
in the past due to relatively low instrument cadences. We report here evidence of such waves directly 
imaged in EUV by the new SDO AIA instrument. In the 2010 August 1 C3.2 flare/CME event, 
we find arc-shaped wave trains of 1-5% intensity variations (lifetime —200 s) that emanate near the 
flare kernel and propagate outward up to —400 Mm along a funnel of coronal loops. Sinusoidal fits 
to a typical wave train indicate a phase velocity of 2200 ± 130kms~^. Similar waves propagating 
in opposite directions are observed in closed loops between two flare ribbons. In the k-uj diagram 
of the Fourier wave power, we find a bright ridge that represents the dispersion relation and can be 
well fitted with a straight line passing through the origin. This k-uo ridge shows a broad frequency 
distribution with indicative power at 5.5, 14.5, and 25.1 mHz. The strongest signal at 5.5 mHz (period 
181 s) temporally coincides with quasi-periodic pulsations of the flare, suggesting a common origin. 
The instantaneous wave energy flux of (0.1-2.6) x 10^ ergs cm~^ s~^ estimated at the coronal base is 
comparable to the steady-state heating requirement of active region loops. 

Subject headings: Sun: activity — Sun: corona — Sun: coronal mass ejections — Sun: flares — Sun: 
oscillations — waves 



1. INTRODUCTION 

In the last decade, observations from SOHO^ TRACE ^ 
Hinode^ and ground-based instruments have led to 
detection of various modes of magnetohydrodynamic 
(MHD) waves in the solar corona (see review by 
iNakariakov fc Verwichtd l2005f ). including (1) oscilla- 
tions or standing waves of slow modes (|Wang et al.l 
120021 : lOfman fc Wan^ 120021). fast kin k modes (peri- 
ods: 2-10 min; lAschwanden et al.l 119991 : iSchriiver et al.l 
119991). and fast sausage modes (periods: 1-60 s; 
INakariakov et all 120031). and (2) propagating waves 
of slow modes (|Of man et al.l 119971: (Deforest & Gur mani 
119981 : iDe MoorteT et al. 2000; Ofman & Wang 2008 ) and 
Alfven waves ([Tomczvk et al.i i2007i: iDe Pontieu et al 
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2007; Ci rtain et al.ll2007l : IOkamoto et al.ll2007l : lJess et al 



2009l : iLiu et a l. 2009; som e of which were int erpreted as 
kink waves, see Van Doo rsselaere et al1l2008[ ). 

Quasi-periodic propagating fast mode magnetosonic 
waves with phase speeds v^\^ — 1000 km s~^ in active re- 
gions remain the least observed among all coronal MHD 
wave s, while single-pulse "EIT waves" ([Thompson et al.l 
119981 ) of typical speeds — 200kms~^ were inte r preted 
as their quiet Sun cou nterpart s (IWu et al.l 120011 : 
Ofman Thompson^ l2002l : cf., IChen Wd [20T¥ . 

first imaged during an eclipse 
2100 kms~^ in a closed loop. 
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iVerwichte eFall (j2005D later observed with TRACE fast 
kink modes of v^\^ = 200-700 km s~^ in an open-field 
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supra- arcade. 

The scarcity of fast wave observations was mainly due 
to instrumental limi tations. The new A tmospheric Imag- 
ing Assembly (AIA: [Lemen et alll2011[ ) on the Solar Dy- 
namics Observatory {SDO) has high cadences up to 10 s, 
short exposures of 0.1-2 s, and a 41' x 41' full- Sun field 
of view (FOV) at V!b resolution, which are all crucial 
for detecting fast propagating features. Within the first 
year of its launch, AIA has detected 10 quasi-periodic 
fast propagati ng (QFP) waves, a mong which the first was 
mentioned by iLiu et all (|2010bl ) and the best example is 
presented here. 

2. OBSERVATIONS AND DATA ANALYSIS 

On 2010 August 1 an eruption (jLiu et al.l l2010al : 
ISchriiver fc Titld 120111 ) occurred in NO A A active region 
11092, involving a coronal mass ejection (CME) and a 
GOES C3.2 flare that started at 07:25 UT and peaked 
at 08:57 UT. 

2.1. Space-time Analysis 
2.1.1. Waves in the Funnel 

In AIA 171 A running difference images (Figure [T]^d) - 
(f). Animation 1(D)) and even direct and base difference 
images (Animations 1(A) and 1(C)), we discovered arc- 
shaped wave trains emanating near the brightest flare 
kernel (box 1 in Figure mb)) and rapidly propagating 
outward along a funnel of coronal loops that subtend an 
angle of '-«^60° near the corona base. They are succes- 
sive, alternating intensity variations of 1-5%, repeatedly 
launched in the wake of the CME during the rise phase 
of the flare (07:45-08:45 UT). The wave fronts continu- 
ously travel beyond the hmb, suggesti ng that the y are not 
propagating over the solar surface like lMoretonI J^GO) or 
EIT waves. They are not observed in the other AIA EUV 
channels, indicating subtle temperature dependence. 
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Fig. 1. — SDO AIA images of QFP waves, (a) 171 A image (see Animation 1(A)) showing the funnel and loop in which fast waves 
propagate, (b) 1600 A image (Animation 1(B)) showing flare ribbons, (c) 171 A base difference image (Animation 1(C)) showing 
dimming behind the CME front. The four brackets mark the smaller FOV of the other panels, (d)-(f) 171 A running difference images 
(Animation l(D))showing successive wave fronts propagating in the funnel. The three curv ed cuts are used to obtain space-time diagrams 
shown in Figure [2] The square box marks the region for Fourier analysis in Section [2.21 (g)-(i) Images of (d)-(f) in the boxed region 
Fourier flltered with a narrow Gaussian centered at the peak in Figure 0Jb) at frequency u = 14.5 mHz (P = 69 s) and wave number 
k = 9.0 X 10~^ Mm"-*^ (A = 110 Mm), which highlight the corresponding QFP wave trains (see Animation 1(E) and Section [2.2. 3|) . 



To analyze wave kinematics, we placed three 20'' 
(14.7 Mm) wide curved cuts that start from the bright- 
est flare kernel and follow the shape of the funnel (Fig- 
ure [TJd)). By averaging pixels across each cut, we ob- 
tained image profiles along it and stacking these profiles 
over time gives space-time diagrams as shown in Figure [21 
where we see two types of moving features: 

(1) The shallow, gradually accelerating stripes repre- 
sent the expanding coronal loops in the CME that have fi- 
nal velocities up to >260 kms~^ as indicated by parabolic 
fits (dashed lines in Figure [2fb)). EUV dimming is evi- 
dent behind these loops (Figures [TJc) and[2fd)), indicat- 
ing evacuation of coronal mass. 

(2) The steep, recurrent stripes result from the arc- 
shaped wave fronts. Sinusoidal fits (Figure [2fe)) to the 
spatial profiles along the central cut yield a projected 
wavelength A = 133 ± 17 Mm and phase velocity = 
2200 ± 130kms~^, giving a period of P = A/vph = 60 ± 
8 s. Linear fits to the space-time stripes from the three 
cuts produced by the same wave front indicate similar 
velocities (Figure [2ja)-(c)). (Such velocities measured 
from projection on the sky plane are lower limits of their 
3D values.) Each wave front travels up to ~400 Mm with 
a lifetime of ~200 s before reaching the edge of AIA's 



FOV, likely resulting from damping and amplitude decay 
with distance (cx 1/r^). 

2.1.2. Waves in Closed Loops 

At the same time, we noticed similar fast propagating 
waves along closed loops between two flare ribbons (Fig- 
ures [TJa) and (b)). The space-time diagram (Figure [5l^a)) 
from the loop-shaped cut reveals steep stripes of both 
positive and negative slopes, particularly near the two 
footpoints, which represent waves propagating in oppo- 
site directions. The bi-directional propagation can be ev- 
idently seen separately in Fourier filtered space-time di- 
agram s (Figures [3fb) and (c); see iTomczvk fc Mclntoshl 
l2009f ). The linearly fitted phase velocities are similar 
in the two directions (1000-2000 kms"^). The sudden 
switches of direction at the western footpoint (top edge 
of the plot) near 08:10 and 08:25 UT suggest wave reflec- 
tion, but a general trend cannot be established. It is thus 
not clear whether the bi-directional waves are generated 
independently, or they are the same wave trains reflected 
repeatedly between the footpoints. 

We find no temporal correlation between the waves in 
the closed loops and those in the funnel that is dominated 
by outgoing waves, except for marginal incoming wave 
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Fig. 2. — Space-time analysis of waves in the funnel, (a)-(c) Running difference space-time diagrams obtained from AIA 171 A images 
along the three cuts shown in Figure md). The insert in (a) offers an enlarged view for the selected region, overlaid with a distance 
averaged profile showing a 43 s periodicity, (d) Base ratio space-time diagram of cut 2 obtained by normalizing image profiles with a 
pre-event profile. All the space-time diagrams are smoothed with a 3 x 3 pixel boxcar, except for (a) which is smoothed in space only, (e) 
Vertical slices of (b) at times and distances marked by the two plus signs. They are snapshots of intensity running difference (x-axis) as 
a function of distance (y-axis) at five consecutive times. Each curve (and thus its average position, marked by the vertical broken line) is 
incrementally shifted by 12 DN which equals AIA's 12 s cadence and thus the x-axis also serves as elapsed time. Each profile is fitted with 
a sine function Asin[27r(r — ro)/A] shown in red, where A is the amplitude, A wavelength, and ro the initial phase in distance. The average 
fitted parameters and their standard deviations are listed. The filled circles mark the delayed occurrences at the average position, to which 
a linear fit indicates a phase velocity -Uph = 2200 ± 130kms~^. (f) Horizontal slices of (d) in the selected region, i.e., temporal profiles of 
intensity base ratio at locations marked by the cross signs. Successive curves at greater distances are shifted upward. The two prominent 
wave periods of 69 and 181 s are marked with slanted lines, indicating wave propagation. 



signals near its base (Figure [2]). Because of their simplic- 
ity (no superposition of bi-directional propagation), we 
choose to further analyze the waves in the funnel with 
Fourier transform as presented below. 



ties of 1630 ± 760kms~^, which cannot be distinguished 
in the observed range up to the Nyquist frequency of 
41.7 mHz given by AIA's 12 s cadence due to the large 
uncertainty. 



2.2. Fourier Analysis of Waves in the Funnel 
2.2.1. Overall Diagram 

We extracted a 3D data cube in (x, t) coordinates, 
i.e., a time series of 171 A running difference images 
for the FOV of Figure HJa) during 07:45-08:45 UT. 
We obtained the Fourier power of the data cube on 
the ij^x") 1 

u) basis of wave number kx and ky and 
frequency v. We then summed the power in the az- 
imuthal direction of cylindrical coordinates {k^O^v)^ 

where k = ^kl + kl fe.g.. iDeForestI 120041 ). This yields 

a k-uo diagram of wave power at a resolution of A/c = 
2.09 X 10"^Mm~^ and = 0.277mHz as shown in 
Figure Hl^a). We find a steep, narrow ridge that de- 
scribes the dispersion relation of the fast propagating 
waves, together with a shallow, diffuse ridge that repre- 
sents those slowly expanding loops at velocities on the 
order of 50kms~^. 

To isolate the fast propagating waves (at the expense of 
reduced frequency resolution), we repeated this analysis 
for a smaller boxed region as shown in Figure [Hd) and a 
shorter duration of 07:58-08:23 UT in which these waves 
are prominent. The resulting k-uj diagram (Figure HJd) 
better shows the steep ridge that can be fitted with a 
straight line passing through the origin. This gives aver- 
age phase ('Uph = v/k) and group {v^^ = du/dk) veloci- 



2.2.2. Temporal Evolution of k^ Diagram 

We repeated the above procedure for a data cube of 
the boxed region during 07:45-08:45 UT masked with a 
running time window that has a full width half maximum 
(FWHM) of 10 minutes with cosine bell tapering on both 
sides. We shifted the window by 1 min at a time (only 6 
such windows are independent in the 1 hr duration) and 
obtained a corresponding k-u diagram, as shown in Fig- 
ures |4l^c)-(f) and Animation 4. The early k-uj diagrams 
are dominated by a shallow ridge with an increasing slope 
that indicates the CME acceleration. When the CME 
front moves out of the FOV, a steep ridge correspond- 
ing to the fast propagating waves becomes progressively 
evident with a slope varying in the 1000-2000 kms~^ 
range. 

2.2.3. Frequency Distribution of Fourier Power 

We note that running difference (time derivative) in 
images used above, similar to a highpass filter, essen- 
tially scales the original signal with frequency u and ap- 
plies a u'^ factor to the Fourier power. To recover the in- 
trinsic power amplitude, we replaced running difference 
images with detrended images obtained by subtracting 
images running-smoothed in time with a 200 s boxcar, 
introducing a low-frequency cutoff of 5 mHz that is be- 
low all strong peaks on the ridge in Figure HJ^b). We 
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Fig. 3. — (a) AIA 171 A running ratio space-time diagram from 
the curved cut along coronal loops shown in Figure [Ha). Distance 
is measured from the eastern footpoint. Note fast waves propa- 
gating in both counter-clockwise (CCW) and clockwise (CW) di- 
rections along the loop, (b) and (c) Fourier filtered version of (a) 
showing CCW and CW waves separately. The linear fits here are 
repeated in (a). 



then repeated the above analysis in Sections 12.2.11 and 
I2.2.2[ The new k-uj diagrams (e.g., Figure Hl^g) vs. (e)) 
exhibit the expected general trend of decreasing power 
with frequency, and as a result the steep ridge becomes 
less evident at high frequencies. 

We averaged the new version (not shown) of the overall 
k-uo diagram of Figure HJ^b) in wave number and obtained 
a power spectrum for the QFP waves (Figure Hl^h)). We 
repeated this for the new k-uo diagrams at different times 
(e.g.. Figure m^g)) and compiled a running spectrogram 
(Figure m^i)). The waves disp l ay a broad frequency dis- 
tribution (cf., iTomczvk et al.l [2QQ7I ). with power peaks 
of ratio 1: 1/4.6 : 1/15.8 (oc zy-i-8±o.2) frequencies 

= 5.5 ± 0.4, 14.5 ± 0.7, and 25.1 ± 0.7 mHz of ratio 
1 : (2.6 ± 0.2) : (4.5 ± 0.3). For comparison, a trian- 
gle wave of vq =5.5 mHz has non-zero Fourier power 
(blue asterisks. Figure HJ^h)) at frequencies of similar ra- 
tio 1 : 3 : 5 that drops faster with 

The Fourier power from running difference and de- 
trended images yields consistent peak frequencies, which 
can be visually identified in the space-time domain. The 
lowest frequency z/q = 5.5 mHz (Pq = 181±13s~3 min) 
manifests as slow modulations in Figures [2fb) and (d) at 
08:06-08:18 UT. The next period 69 ± 3s (14.5 mHz), 
dominating the power from running difference images 



(Figure HJ^b)), matches the temporal spacing between 
bright stripes near 08:08 UT in Figures [2f a)- (c) and the 
period given by the sinusoidal fits (Figure [2fe)). The 
corresponding wave fronts are prominent in the origi- 
nal and Fourier filtered images (Figures [Hd)-(i), Ani- 
mation 1(E)). These two periods are also evident in the 
emission profiles of Figure EJf). The higher frequency 
25.1 mHz (40 ± 1 s) has considerably weaker power and 
a close frequency of 23 mHz (see Figure lU^d)) can be 
seen in the spacing of narrow stripes near 08:01 UT (Fig- 
ure Efa)), when the other two frequencies are not yet 
strong. 

2.3. Common 3 min Periodicity in Waves and Flare 
As shown in Figure [5jb), the RHESSI X.-Tdiy flux and 
AIA 1600 A fluxes of flare ribbons (particularly the 
brightest one in box 1 where the funnel is rooted, see 
Figure [TJb) and Animation 1(B)) exhibit bursty bumps 
at a 3 min period (5.5 mHz). The onsets of these pul- 
sations (vertical dotted lines) coincide with those of the 
slow modulations on the QFP waves (Figure [5fa)). This 
can be also seen in the wavelet power of these flare emis- 
sions (Figures [5fd)-(g)). The Fourier power of the X-ray 
flux (green curve. Figure lU^h)) is consistent with that 
of the QFP waves <10 mHz, but significantly lower at 
higher frequencies. It also matches that of the triangle 
wave up to the third harmonics because of its triangular 
pulse shape (Figure [5fc)). In contrast, the 1600 A flux of 
a background plage (box 3 in Figure [TJb)) is constantly 
dominated by the 5 min (3.5 mHz) photospheric p-mode 
oscillations (Figure [5jb) and (f)). 

2.4. Estimate of Wave Energy and Magnetic Field 

The energy flux carried by the QFP waves can be esti- 
mated with the kinetic energy of th e perturbed plasma , 
E = p{5vfv^i,/2 > p{5I/I)^vlJ8 (|Aschwandenll200l . 
where we have assumed that the observed intensity vari- 
ation 61 results from density modulation 6p and used 
Sv/vph ^ ^p/p = ^^/(2/) for magnetosonic waves since 
/ oc p^. If we take v^^ = 1600 kms~^ and SI /I = 1%- 
5% observed in the mid-range of the funnel (200 Mm 
from the flare kernel), and use the corresponding num- 
ber density Ue ^ 10^ cm~^ estimated with the 17 1 A 
channel response (following iDe Pontieu et al]l2011f ). we 
reach an energy flux E > (0.1-2.6) x 10^ ergscm"^ s~^. 
The diameter of the funnel here has increased ~10 times 
from the coronal base, where the the wave energy flux 
shall be ^10^ times higher by continuity of energy flow, 
if we assume the waves being generated there and con- 
sider damping on their path. This energy flux is more 
than sufficient for heatin g the local active region loops 
(jWithbroe fc Novesll 19771 ). However, considering the lim- 
ited temporal and spatial extent of these waves, they are 
unlikely to play an important role in heating the quies- 
cent global corona. 

Assuming the measured phase speed equal to the 
fast mode speed along magnetic field lines in the funnel, 
which is the Alfven speed va = B/^Airp, the magnetic 
field strength is estimated as P = Vph^/47^p ^ 8 G. 

3. DISCUSSION 

We propose that these QFP waves imaged with AIA's 
unprecedented capabilities are fast mode magnetosonic 



Fast Mode MHD Waves Observed by SDO AIA 



5 



, (o) powei;: fujl,FOV, ,07:4^-08:^5 



(b),Powejr: In ,box, ,07:5^-08:^5 (h) Powj^r spect^-um (i), Spec^rogrqrp (dynprnic^pec,) 

^3-pomt boxcorJ 



2 4 6 8 
Wove number k (10"' 



: Linear fit to ridge: 
; Vp,= 1630±760 km/s 
: <afc>=0.63x10"^ Mm" 



: Two power peaks: 



ife= (0.9 ±0.46)x 10-2 Mm- 
X=110±60 Mm 



ib= (0.45 ±0.57)x 10-2 Mm- 
X=220±2ff0 Mm 



RHESSI X-royl 




25.1 ±0.7 mHz 
4d±1 s 



14.5±0.7 mHz 
69 ± 3 s 



9 
o 

^-4.0 



10 
Mm"') 



2 4 6 8 
Wove number k (10"' 



10 
Mm"') 



0.010 0.100 
Power (ON^) 



07:50 08:00 08:10 08:20 08:30 
Time (UT, since 2010/08/01 07:47:36) 



(c) 07:45:30-07:55:30 , (d,) 07:55:30-08:05:30 , (e) 08:05:30-08:15:30 



(f) 08:15:30-08:25:30 (g) Log e) by detrended ,doto 



^30- 
E 



/ O 25.1 ±0.7 mHz ; § 

/ o / o 

I 14.5±0.7 mHz ' 



5.5 ±0.4 mHz i 




02468 10 12 02468 10 12 02468 10 12 02468 10 12 02468 10 12 
Wove number k (10"' Mm"') 

Fig. 4. — Fourier analysis of wave properties, (a) Fourier power {k-uo diagram) of a 3D data cube of 171 A running difference images 
during 07:45-08:45 UT in the full FOV of Figure [fta). (b) Same as (a) but for 07:58-08:23 UT in the boxed region of Figure [Hd) . The 
white dashed line (repeated in (a)) is a power- weighted linear fit to the (fc, u) positions of pixels greater than 10% of the maximum power 
in the k < /cmax range (marked by the vertical dotted line), (c)-(f) Same as (b) but for images masked with a running time window whose 
FWHM is labeled on eac h panel (see Animation 4). (g) Same as (e) but on log scale from detrended (rather than running difference) 
images (see Section r2.2.3|) . The diffuse horizontal band is an artifact at the 5 mHz detrending cutoff frequency and is ~5% of the QFP 
wave power here, (h) Power spectrum vs. frequency obtained by averaging in k < /cmax on a k-uu diagram that is the same as (b) but from 
detrended images, (i) Spectrogram obtained by compiling wave number averaged power as shown in (h) from k-uu diagrams at different 
times as shown in (g). The x-axis here refers to central times of the running window. Prominent "islands" are contoured at the 50% level; 
their peaks are marked by plus signs and the peak frequencies (periods) by horizontal dotted lines. The frequency uncertainties are the 
standard deviations within the contours. 



waves that have been theoretically predicted and sim- 
ulated fe.g.. iBogdan et al.l[2QQl iFedun et aIll2Qlll : Of- 
man et al. in preparation), but rarely observationally de- 
tected. We speculate their possible origin as follows. 



can modulate ener gy release and flare emission 
(|Foullon et al.l [20051 ). This periodicity is the same 
as that of 3 min chromospheric oscillations, fur- 



1 



The accompanying CME is unlikely to be the 
wave trigger, because it takes place gradually for 
~30 min wave periods. Figure El^b)) and its 
single pulse would be difficult to sustain oscilla- 
tions lasting hr as observed here without being 
damped. However, the environment in its wake 
might be favorable for these waves. 



nection JChen & Priest 


120061; iHessland et al.||2009|; 


iMcLauehlin et al.ll2009 


)• 



2. The common 3 min periodicity (Section 
of the QFP wa ves and flare quasi-periodic 
pulsations (QPPs: iNakariakov fc Melnikovl 120091 : 



ikuprivanova et al.l l2010f ) suggests a common ori- 
gin. Quasi-periodic magnetic reconnection and 
energy release can excite both flare pulsa tions 
(jOfman fc Suil 120061 : [Fleishman et al.l I2008D and 



MHD oscillations that drive QFP waves, or in 
turn, MHD oscillations responsible for the waves 



However, the deficit of flare power at higher wave 
frequencies (^10 mHz, Figure [D^h)) is somewhat 
puzzling. Perhaps the waves are driven by a multi- 
periodic exciter that produces no detectable flare 
signals at these frequencies. A future study of sim- 
ilar events will further shed light on the nature of 
these waves. 
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